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Abstract 

Although the MiniBooNE experiment has severely restricted the possible ex- 
istence of light sterile neutrinos, a few anomalies persist in oscillation data, and 
the possibility of extra light species contributing as a subdominant hot (or warm) 
component is still interesting. In many models, this species would be in thermal 
equilibrium in the early universe and share the same temperature as active neutri- 
nos, but this is not necessarily the case. In this work, we fit up-to-date cosmological 
data with an extended ACDM model, including light relics with a mass typically in 
the range 0.1-10 eV. We provide, first, some nearly model-independent constraints 
on their current density and velocity dispersion, and second, some constraints on 
their mass, assuming that they consist either in early decoupled thermal relics, or in 
non-resonantly produced sterile neutrinos. Our results can be used for constraining 
most particle-physics-motivated models with three active neutrinos and one extra 
light species. For instance, we find that at the So" confidence level, a sterile neutrino 
with mass nris = 2 eV can be accommodated with the data provided that it is ther- 
mally distributed with Tg/T^ < 0.8, or non-resonantly produced with AAofj < 0.5. 
The bounds become dramatically tighter when the mass increases. For < 0.9 
eV and at the same confidence level, the data is still compatible with a standard 
thermalized neutrino. 



1 Introduction 

Neutrino oscillation is a well studied phenomenon, confirmed by strong experimental 
evidences. Most experimental results are well explained with a three-neutrino oscillation 
model, involving two independent and well-measured square-mass differences: Am^^j = 
(7.59 ±0.21) X 10-5 eV^ [1] and /S.matm = (2.741°;^^) x lO'^ eV^ [2]. However, some other 
experiments have shown some anomalies which do not fit in this hypothesis (LSND [3], 
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Gallium experiments [Ij, MiniBooNE low energy anomaly [5]). These anomalous results 
might be due to unknown systematic effects, but all attempts to identify such systematics 
have failed until now. Otherwise, they could be interpreted as exotic neutrino physics. 

In Ref. j6j, the MiniBooNE anomaly was explained through a renormalization of the 
absolute neutrino flux and a simultaneous disappearance of electron neutrinos oscillat- 
ing into sterile neutrinos (with Pu^^ue = ^-^^-om)- The LSND and Gallium radioactive 
source experiment [TJ El [9] anomalies have been studied in Ref. [TO] , where is it claimed 
that all these anomalies could be interpreted as an indication of the presence of, at least, 
one sterile neutrino with rather large mass (few eV's). Ref. [11] also studied the com- 
patibility of the Gallium results with the Bugey and Chooz [13] reactor experimental 
data, concluding that such a sterile neutrino should have a mass between one and two 
eV's. Finally, the MiniBooNE collaboration performed global fits of MiniBooNE, LSND, 
KARMEN2, and Bugey experiments in presence of a fourth sterile neutrino [H] (assuming 
no renormalization issue for MiniBooNe unlike Ref. [6]). When all four experiments are 
combined, the compatibility between them is found to be very low (4%); however, when 
only three of them are included, the compatibility level is usually reasonable (the largest 
tension being found between LSND and Bugey). In this analysis, the preferred value of 
the sterile neutrino is usually smaller than leV, but still of possible cosmological relevance 
(for instance, for all four experiments, the best fit corresponds to Am^ ~ 0.2 — 0.3 eV^). 

These various developments suggest that it is important to scrutinize cosmological 
bounds on scenarios with one light sterile neutrino, which could help ruling them out, 
given that current bounds on the total neutrino mass assuming just three active neutrinos 
are as low as ^2^^ < 0.61eV (using WMAP5, BAO and SN data [15]). This result 
cannot be readily applied to the models which we consider here. Indeed, scenarios with 
extra neutrinos require a specific cosmological analysis, for the simple reason that besides 
affecting the total neutrino mass, additional neutrinos also increase the abundance of 
relativistic particles in the early universe. 

From the point of view of Cosmology, there have been many works constraining si- 
multaneously the sum of neutrino masses and the contribution to the relativistic energy 
density component of the Universe, parametrized as the effective number of neutrinos, iVefj 
(see for example [TUlIITlIIHlIinilSnilSJ)- Most of these works assume either that the heav- 
iest neutrino (and hence the most relevant one from the point of view of free-streaming) 
has a thermal distribution, sharing the same value of temperature as active neutrinos, or 
that all neutrinos are degenerate in mass. However, the results of Refs. [221123] can also be 
applied to the case of very light active neutrinos plus one heavier, non-necessarily thermal 
sterile neutrino, which is the most interesting case for explaining oscillation anomalies. In 
terms of physical motivations, it is very likely that the light sterile neutrino required by the 
LSND anomaly acquires a thermal distribution in the early universe, through oscillations 
with active neutrinos in presence of a large mixing angle [2^. On the other hand, there 
are some proposals to avoid these contrains (for a list of some scenarios, see [25]). One of 
such possibilities is based on a low reheating temperature (Tr) Universe [2Sl EZl EHl 1^ - 
in which, for a sufficiently low Tr, the sterile neutrinos could be non-thermal [SO] and its 
production would be suppressed [28] , such that usual cosmological bounds are evaded. In 
fact, in these models, sterile neutrinos are allowed to have a large mass without entering 
in conflict with other experimental results, while T^ec ^ ^j? ^ lOMeV (Tdec being the 
temperature of the cosmic plasma at neutrino decoupling). 
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In absence of thermalization, cosmological bounds on the sterile neutrino mass become 
potentially weaker. Hence, it is interesting to study the compatibility of recently proposed 
scenarios with a light sterile neutrino with the most recent cosmological data, keeping 
in mind the possibility of a non-thermal distribution. The goal of this paper is hence 
to study the compatibility of cosmological experimental data (WMAP5 plus small-scale 
CMB data, SDSS LRG data, SNIa data from SNLS and conservative Lya data from VHS) 
with the hypothesis of a sterile neutrino with the characteristics sketched above, i.e., with 
a mass roughly of the order of the electron- Volt, and a contribution to Ncs smaller than 
one. 



2 Light sterile neutrino in cosmology: physical ef- 
fects and parametrization 

If a population of free-streaming particles becomes non-relativistic after photon decou- 
pling, its physical effects on the cosmological background and perturbation evolution are 
mainly described by three quantities: 

1. its contribution to the relativistic density before photon decoupling, which affects 
the redshift of radiation/matter equality, usually parametrized by an effective neu- 
trino number (standing for the relativistic density of the species divided by that of 
one massless neutrino family in the instantaneous decoupling (id) limit): 



^rel 



dp P^fip) 



I 



815 



(1) 



with Tf = {4/liy/%, 



its current energy density, which affects (i) the current energy budget of the Universe 
(with various consequences for the CMB and LSS spectra, depending on which 
other parameters are kept fixed), and (ii) the amplitude reduction in the small- 
scale matter power spectrum due to these extra massive free-streaming particles, 
parametrized by the dimensionless number utg: 



^ Idpp'fip) 



(2) 



where p° is the critical density today and h the reduced Hubble parameter. 



3. the comoving free-streaming length of these particles when they become non-relativistic, 
which controls the scale at which the suppression of small-scale matter fluctuations 
occurs. This length can easily be related to the average velocity of the particles 

today, {vs 



^The minimum comoving free-streaming wavenumber kfs is controlled by f2„i and by the ratio 
fi(^nr)/(ws(^nr)) evaluated when T = m, i.e. when a(inr) ~ {vs{to))o-{ta)- Given that (ws(inr)) ^ 
{vs{to))a{to) / a{tnr) , the minimum comoving free-streaming length just depends on {vs(to)) and flm- 
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However, for whatever assumption concerning the phase-space distribution function f{p), 
the three numbers (AiVefj, Ug, (fs)) satisfy a constraint equation. Indeed, the average 
velocity of the particles today (assumed to be in the non-relativistic regime) is given 
exactly by 

(vs) = r o , ;/ = — — = 5.618 X 10~^ — 3 

in units where c = ks = h = 1, and taking Tcmb = 2.726K. Hence, the three physical 
effects described above depend on only two independent parameters. 

Reducing the physical impact of any population of massive free-streaming particles to 
these three effects (and two independent parameters) is a simplification: two models based 
on different non-thermal phase-space distributions f{p) can in principle share the same 
numbers {ANes, i^s, (vs)) and impact the matter power spectrum differently. Indeed, the 
free-streaming effect depends on the details of f{p) (including high statistical momenta 
like Jdpp^f{p), etc.) However, the conclusions of Ref. [TH] indicate that for many 
models with non-thermal distorsions, observable effects can indeed be parametrized by 
two combinations of (AA^'eff, ujg, {v^)) with good accuracy: other independent parameters 
would be very difficult to observe!. 

Let us compute the three parameters (AA^^efr, ujg, {vs)) for simple cases. For one species 
of thermalized free-streaming particles with mass m^, sharing the same temperature as 
active neutrinos in the instantaneous decoupling limit, one gets: 

TUs , , 771^ T'}f 0.5283 meV 

94.05 eV 180C(3) m<j 

For a light thermal relic with a Fermi-Dirac distribution and a different temperature Tg, 
these quantities become 

AM f^A' fTsY 0.5283 meV fT,\ 

For a non-thermal relic with a free function f{p), there is an infinity of possible models. 
A popular one is the Dodelson-Widrow scenario [32] (also referred as the "non-resonant 
production scenario"), motivated by early active-sterile neutrino oscillations in the limit 
of small mixing angle and zero leptonic asymmetry, which corresponds to the phase-space 
distribution 

where x is an arbitrary normalization factor. In this case, in the approximation = T^f, 
the three "observable" parameters read 

ms , , 0.5283 meV 

94.05eV m. 



^This conclusion does not apply when the non-thermal distribution f{p) has a sharp peak close to 
p = 0. In this case, particles with very small momentum should be counted within the CDM component, 
not within the extra massive free-streaming component. Otherwise, one would obtain values of LUg and 
(vs) based on an averages between cold and hot /warm particles; then, these parameters would not 
capture the correct physical effects (see [3T]) 
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Hence, a Dodelson-Widrow (DW) model shares that same "observable" parameters (AA^eS, 
cUs, (vs)) as a thermal model with TTjthermai _ ^dw^i/a ^^^^ j,^ _ ^1/4 ji^^ Actually, for 

these two models, the degeneracy is exact: it can be shown by a change of variable in 
the background and linear perturbation equations that the two models are strictly equiv- 
alent from the point of view of cosmological observables [331 [TH] . As mentioned before, 
in the general case, two models sharing the same (AA'eg, Ug, (vs)) are not always strictly 
equivalent, but can be thought to be hardly distinguishable even with future cosmological 
data. For instance, the low-temperature reheating model analyzed in [271 [28] leads to a 
distribution of the form 

This model would in principle deserve a specific analysis, but in good approximation we 
can expect that by only exploring the parameter space of thermal models (or equivalently, 
of DW models), we will obtain some very generic results, covering in good approximation 
most possibilities for the non-thermal distorsions. 



3 Analysis 

3.1 Data 

In the following sections, we will present the results of various runs based on the Boltz- 
mann code CAMB [HI] and cosmological parameter extraction code CosmoMC [35]. We 
modified CAMB in order to implement the proper phase-space distribution f{p) of the 
thermal or DW model. For simplicity, we assumed in all runs that the three active neu- 
trinos can be described as massless particles. In order to obtain a Bayesian probability 
distribution for each cosmological parameters, we ran CosmoMC with flat priors on the 
usual set of six parameters Ub, cudm = ujg + i^cdm, r, Ag, Ug (see e.g. [36]), plus two 
extra parameters describing the sterile neutrino sector, that will be described in the next 
sections. We choose the following data set: WMAP5 [37] plus small-scale CMB data 
(ACBAR [38], CBI [39], Boomerang [50]), the galaxy power spectrum of the SDSS LRG 
|41j with flat prior on Q [121 S3j, SNIa data from SNLS [H] and conservative Lyman-a 
data from VHS [45j . We do not include more recent Lyman-a data sets, which have much 
smaller errorbars, but for which the deconvolution of non-linear effects depends on each 
particular cosmological model, and requires specific hydrodynamical simulations. 

3.2 General analysis 

Our first goal is to obtain simple results with a wide range of applications. Hence, we 
should not parametrize the effect of sterile neutrinos with e.g. their mass or temperature: 
in that case, our results would strongly depend on underlying assumptions for f{p). It 
is clear from section [2] that nearly "universal" results can be obtained by employing 
two combinations of the "observable parameters" AA^^eff, oJs and (f^) (and eventually of 
other parameters of the ACDM model). Here we choose to vary the current dark matter 
density fraction fs = Us/ioJs + t^cdm) and the current velocity dispersion {vs). As will be 
clear from our results, these two parameters capture the dominant observable effects, and 
lead to very clear bounds, since their correlation with other ACDM model parameters 
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is insignificant. Our limits on fs and (vs) apply exactly to the thermal case and DW 
case, and approximately to most other cases (modulo the caveat described in the second 
footnote of section El). 

Our parameter space is represented in Figure [TJ We adopt a logarithmic scale for (vs) 
and display the interesting range 

1 km/s < (vs) < 1000 km/s . (9) 

Indeed, with out dataset, particles with smaller velocities would be indistinguishable from 
cold dark matter; instead, particles with larger velocities would either have AN^s > 1 (a 
case beyond the motivations of this work, and anyway very constrained by the data) or 
fs < 0.02 (being indistinguishable from extra relativistic degrees of freedom). Assuming 
a particular value for u^m = + i^cdm and for ANes, it is possible to compute the 
velocity dispersion (vs) as a function of fs- Since the CMB and LSS data give precise 
constraints on Udm, regions of equal ANea correspond to thin bands in the {fs, (vs)) 
plane. We show these bands in Figure [T] for 10"'^ < AiVcfj < 1 under the assumption that 
^dia = 0.11 ± 0.01, which corresponds roughly to the 95% confidence limits (CL.) from 
all our runs. These iso-AA^^efr bands are completely model-independent. 

Instead, regions of equal mass can only be plotted for a particular model. In Figure [H, 
we show the bands corresponding to m =1 eV and 10 eV, either in the case of early 
decoupled thermal relics (blue/dotted hues) or in the DW case (green/dashed lines). For 
any given mass, these bands intersect each other in a location corresponding to the case 
of one fourth standard neutrino species with AA^eS = 1- 

We ran CosmoMC with top-hat priors on fs (in the physical range [0,1]) and on 
log]^Q[(fs)/lkm/s] (in the range [0,3] motivated by the previous discussion). Our results 
are summarized in Figure [1] (bottom). We see that the upper bound on fs decreases 
smoothly as the velocity dispersion increases: when the particles have a larger velocity 
dispersion, their free-streaming wavelength is larger, so the step-like suppression in the 
power spectrum (which amplitude depends on fs) is more constrained. For (vs) ~ 1 km/s, 
we find /, < 0.1 at the 2a C.L., while for {vs) ~ 100 km/s, we find /, < 0.06 at the 2a 
CL. When the velocity dispersion becomes larger than 100 km/s, the upper bound on 
fs decreases even faster as a function of (vs). This is the case of a HDM component with 
significant contribution to the number of relativistic d.o.f., for which the observational 
bounds derive from a combination of the first and second effects described in section [2j 
in this limit, in addition to being sensitive to the free-streaming effect, the data disfavor 
a significant increase of the total radiation density corresponding to AA^^cfi of order one 
or larger. 

We should stress that the details of our results depend on the underlying priors. For 
instance, one could use a flat prior on (vs) instead of its logarithm. Running in the range 
< (vs) < 1000 km/s with such a prior would give more focus on the large- (f^) allowed 
region of Figure [H However, it would be more interesting to focus on small velocities, 
in order to understand how our results can be extended without any discontinuity to the 
case of warmer and heavier dark matter. For this purpose, we ran CosmoMC with a 
top-hat prior on < (vs) < 1 km/s, and obtained the results shown in Figure O These 
results are identical to those published in Reference |16] (figure 7). By gluing figure [H 
on top of [2], one can obtain a full coverage of the parameter space of ACDM models 
completed by one extra (hot or warm) dark matter species. Figure [2] shows the transition 
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from the region in which this extra species is indistinguishable from cold dark matter 
(when (vs) < 0.1 km/s, the fraction fg is unconstrained) to the region in which it is 
warm (for 0.4 < (vg) < 1 km/s, there is a nearly constant bound fg ^ 0.1 at the 2-a 
level) . Figure [T] shows instead the transition from warm particles to hot particles (with 
velocities comparable to those of active neutrinos). The two plots perfectly match each 
other along the (vg) = 1 km/s axis, on which the sterile neutrino fraction is bounded by 
/s<0.1 (2-a). 

3.3 Mass/temperature bounds in the thermal case 

We now focus on the particular case of early decoupled thermal relics, with a Fermi- 
Dirac distribution and a temperature Tg. These models can be parametrized by the mass 
rrig and the temperature in units of the neutrino temperature, Tg/T^f. Our parameter 
space - and the correspondence with the previous parameters AN^s, fg, (vg) - is shown 
in Figure [31 In this analysis, we want to focus again on light sterile neutrinos rather 
than WDM; hence we are not interested in velocities smaller than 1 km/s today. We 
are not interested either in the case of enhanced particles with AA^efr > 1- Then, as can 
be checked in Figure [Sj the ensemble of interesting models can be covered by taking a 
top-hat prior on logiQ^rrig/l eV) in the range [—1, 2], and on Tg/T^ in the range [0, 1]. 

The likelihood contours obtained for this case are shown in Figure O (bottom). They 
are consistent with our previous results: when AN^a ~ 10~^ (and hence Tg/T^f ~ 0.3), 
the upper bound on the sterile neutrino fraction is fg < 0.1 at the 2a C.L.; then this bound 
decreases smoothly when Tg increases. For a fourth standard neutrino with Tg = T^'^, the 
2a C.L. (resp. 3a C.L.) bound is < 0.4 eV (resp. 0.9 eV). 

This figure can be conveniently used for model building: for a given value of the 
mass, it shows what should be the maximal temperature of the thermal relics in order to 
cope with cosmological observations; knowing this information and assuming a particular 
extension of the particle physics standard model, one can derive limits on the decoupling 
time of the particle. For instance, for a mass of rUg = 0.5 eV one gets Tg/T^ < 0.9; for 
nig = 1 eV, Tg/T^ < 0.7; while for = 5 eV, Tg/T^ < 0.5. This figure can also be 
applied to thermally produced axions, like in Refs. [171 HH]. 

3.4 Mass bounds in the DW case 

Finally, for Dodelson-Widrow relics with a distribution function equal to that of stan- 
dard neutrinos suppressed by a factor x (which is equal by definition to AA^efi), we can 
parametrize the ensemble of models by rUg and x- parameter space - and the cor- 
respondence with fg, (vg) - is shown in Figure HI Like in the previous section, we are 
not interested in a current velocity dispersion smaller than 1 km/s today. Then, as can 
be checked in Figure HI the ensemble of interesting models can be covered by taking a 
top-hat prior on log^gl^s/l range [—1,2]; in this range, values of x smaller 

than 10^^ would correspond to tiny values of fg, i.e. to particles indistinguishable from 
massless particles; so, we can take a flat prior on logio(x) i^i the range [—2, 0]. 

The likelihood contours obtained for this case are shown in Figure HI (bottom). We 
are not surprised to find once more an allowed region corresponding to < 0.1 at the 2a 
C.L. when AiVefj = x ~ 10~^ is negligible with respect to one, or less when AiVeff grows 
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closer to one. For a fourth standard neutrino with Tg = , the two definitions of the 
mass (following from the thermal or from the DW cases) are equivalent, and indeed we 
find rus < 0.4 eV (2a C.L.) or < 0.9 eV {Sa C.L.) like in section [3Jl 

This figure can also be useful for model building: for a given value of the mass, it 
shows what should be the maximal value of x compatible with cosmological observations; 
in turn, this information can be used to put bounds on the mixing angle between this 
relic and active neutrinos in non-resonant production models d la Dodelson Sz Widrow. 
For instance, for a mass of rris = 1 eV, the 2a C.L. gives x ^ 0.5; for = 2 eV, we get 
X < 0.2; while for = 5 eV, we get x ^ 0.1. 

3.5 Comparison with previous work 

The ensemble of cosmological models that we are exploring here is not different from that 
studied by Dodelson, Melchiorri & Slosar [23] (called later DMS) or by Cirelli & Strumia 
[^2] (called later CS); the difference between these works and the present analysis consists 
in a different choice of parameters, priors, data set, and also methodology in the case of 
CS. 

For instance. Figure 6a of CS presents constrains in the space (logiQANes, logig^s) 
assuming a DW scenario. Hence, their parameter space is identical to the one we used 
in section excepted for the prior range (which is wider in their case). As far as the 
data set is concerned, CS use some CMB and galaxy spectrum measurements which are 
slightly obsolete by now; on the other hand, they employ some additional information 
derived from BAO experiments, and use SDSS Lyman-a data points that we conserva- 
tively excluded from this analysis, since they assume a ACDM cosmology. Finally, CS 
performed a frequentist analysis, and their bounds are obtained by minimizing the 
over extra parameters (while in the present Bayesian analysis, we marginalize over them 
given the priors). 

In order to compare our results with CS, we performed a run with top-hat priors on 
logio X = logio(^^eff) in the range [—3, 1], and on log^Q^nis/l eV) in the range [—1, 3]. In 
this particular run we compute the 90%, 99% and 99.9% C.L., following CS. Our results 
are shown in Figure El and are consistent with those of our general analysis. 

In spite of the different data set and methodology, the 90% and 99% contours are 
found to be in very good agreement with CS in most of the parameter space. The major 
difference lies in the small mass region, for which CS get more conservative limits on AA^^efr 
than we do, and find a preference for non-zero values of the effective neutrino number 
0.5 < AA^eff < 4 (at the 90% C.L.). This qualitative behavior has been nicely explained 
in Refs. US] . It is due to the non-linear corrections applied to the theoretical linear 
power spectrum before comparing it with the observed SDSS and 2dF galaxy power 
spectra. The approach used in this work (and in the default version of CosmoMC) 
consists in marginalizing over a nuisance parameter Q (describing the scale-dependence 
of the bias) with a flat prior. Instead, following Ref. [IH], CS impose a gaussian prior on 
Q. This results in biasing the results towards larger values of Aefr, and flnding marginal 
evidence for AA'efr > 0. Of course, this assumption might turn out to be correct; however, 
it is argued in Refs. [HI 03] that our knowledge on Q (based essentially on N-body 
simulations for some particular cosmological models) is still too uncertain for getting 
deflnite predictions. 
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The analysis of DMS is Bayesian, like ours. The authors use top-hat priors on the two 
parameters —3 < log^Q^nis/l eV) < 1 and < < I, roughly the same data set as CS, 
and employ the distribution function of early decoupled thermal relics. Our results based 
on the same priors (but a different data set) are shown in Figure [HI and are consistent 
with the previous sections: at the 2a level, Ug is such that fg < 0.1 for AiVefj ~ 10~^; 
then, the bound on fg (and therefore on Ug) decreases smoothly when decreases (and 
therefore (vg) increases). These results differ significantly from those of DMS, who find 
that the upper bound on Ug peaks near m ~ 0.25 eV and then decreases quickly. We 
do not observe such a behavior: our upper bound on Ug increases (not so smoothly, but 
still monotonically) when mg increases, in agreement with all previous results in this 
paper. This difference is most likely due to the use made by DMS of more aggressive 
Lyman-a data from SDSS, of different galaxy power spectrum data, and of a prior on Q, 
as in CS. This data set puts stronger limits on a possible suppression of the small scale 
power spectrum. Actually, in absence of sterile neutrinos, the same combination of data 
is known to produce very strong bounds on neutrino masses, and to prefer AA^gg slightly 
larger than one 03]; in presence of light sterile neutrinos, the results of DMS show that 
this data also imposes a strong bound Ug < 0.001 for leV < rris < 10 eV, due to its 
sensitivity to the sterile neutrino free-streaming effect. Our large scale structure data set 
(conservative Lyman-a data from VHS, SDSS-LRG and flat prior on Q) is not able to 
exclude this region. 

4 Conclusions 

In this work, we studied the compatibility of cosmological experimental data with the 
hypothesis of a non-thermal sterile neutrino with a mass in the range 0.1 — 10 eV (or 
more), and a contribution to iVefj smaller than one. We computed Bayesian confidence 
limits on different sets of parameters, adapted to the case of thermal relics (section [33]), of 
non- resonantly produced sterile neutrinos d la Dodelson & Widrow (DW, section [3^ . or 
of generic parameters leading to nearly model-independent results (section 13.21) . In each 
case, we performed a specific parameter extraction from scratch, in order to obtain reliable 
results assuming flat priors on the displayed parameters. For simplicity, we assumed that 
the masses of the three active neutrinos are negligible with respect to that of the sterile 
neutrino. 

For a cosmological data set consisting in recent CMB and LSS data, as well as older 
but very conservative Lyman-a data, we found the conditional probability e.g. on the 
mass of a thermal relic given its temperature, or on the mass of a DW neutrino given its 
density suppression factor, etc. These proabilities are such that if the fourth neutrino is a 
standard one (with AN^s = 1), it should have a mass nig < 0.4 eV {2a C.L.) or nig < 0.9 
eV {3a C.L.). 

At the 3a C.L., a mass nig = 1 eV can be accommodated with the data provided that 
this neutrino is thermally distributed with Tg/T^ < 0.97, or non-resonantly produced 
with AAefT < 0.9. The bounds become dramatically tighter when the mass increases. 
At the same confidence level, a mass of just ms = 2 eV requires either Tg/T^^ < 0.8 or 
AA'eff < 0.5, while a mass = 5 eV requires Tg/T^ < 0.6 or AAeff ^ 0.2. 

Our bounds can hopefully be used for constraining particle-physics-motivated models 
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with three active and one sterile neutrinos, as those investigated recently in order to 
explain possible anomalies in neutrino oscillation data. Many of these models can be 
immediately localized in our figures [3] or HI For sterile neutrinos or other particles which 
do not fall in the thermal or DW category, a good approximation consists in computing 
their velocity dispersion and localizing the model in our figure Future neutrino oscil- 
lation experiments are expected to test the self-consistency of the standard three-neutrino 
scenario with increasing accuracy. If anomalies and indications for sterile neutrinos tend 
to persist, it will be particularly useful to perform joint analysis of oscillation and cos- 
mo logical data, using the lines of this work for the latter part. 
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Figure 1: (Top) the parameter space {fs,{vs)) chosen in our general analysis. The 
thin bands delimited by red/solid lines show regions of equal AN^s (assuming cjdm = 
O.lliO.Ol); these bands are fully model-independent. We also show the model-dependent 
regions of equal mass, delimited by blue/dotted lines for the case of early decoupled ther- 
mal relics, and consisting in horizontal green/dashed lines for Dodelson-Widrow sterile 
neutrinos. (Bottom) same with, in addition, the regions allowed at the 68.3% (lex), 95.4% 
{2a) and 99.7% {3a) CL. by our cosmological data set, in a Bayesian analysis with fiat 
priors on fs and \ogiQ{vs) within the displayed range. 
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Figure 2: la, 2a and 3cr contours of the marginalized likelihood for the two parameters 
(/s! {'^s))i with different priors than in previous figures. As explained in the text, this plot 
shows the region where the sterile neutrino is heavy and behaves like warm dark matter, 
in complement to Figure 1, which is based on a different range/prior for (f^) adapted to 
the case of a light, hot sterile neutrino. 
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Figure 3: (Top) the parameter space (m^jT^/T^'^) used in the particular case of early 
decoupled thermal relics of temperature (with T^*^ = (4/ll)^/^T^). The thin bands 
delimited by blue/dot-dashed lines show regions of equal fs (assuming cjdm = 0.11 ±0.01); 
the magenta/dotted lines correspond to fixed values of the velocity dispersion today; 
horizontal red/solid lines to fixed AA^^eff. (Bottom) same with, in addition, the regions 
allowed at the 68.3% (lex), 95.4% {2a) and 99.7% (3(t) C.L. by our cosmological data 
set, in a Bayesian analysis with fiat priors on log^gl^^s) Ts/T^ within the displayed 
range. 
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Figure 4: (Top) the parameter space {msjx) used in the particular case of DW rehcs. The 
thin bands dehmited by blue/dot-dashed hues show regions of equal fs (assuming = 
O.lliO.Ol); the magenta/dotted lines correspond to fixed values of the velocity dispersion 
today; horizontal red/solid lines to fixed AA'cff. (Bottom) same with, in addition, the 
regions allowed at the 68.3% (lo"), 95.4% {2a) and 99.7% {3a) C.L. by our cosmological 
data set, in a Bayesian analysis with flat priors on logiQ(ms) and logiQ{x) within the 
displayed range. 
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Figure 5: (To'p) the parameter space (AAi'effjms) used for comparison -with Cirelli & 
Strumia in the particular case of DW rehcs. The thin bands dehmited by blue/dot- 
dashed hues sho'w regions of equal fs (assuming ujAm = 0.11 ± 0.01); the magenta/dotted 
lines correspond to fixed values of the velocity dispersion today; horizontal red/solid lines 
to fixed ANcs- (Bottom) same with, in addition, the regions allowed at the 90%, 99% 
and 99.9% C.L. by our cosmological data set, in a Bayesian analysis with flat priors on 
log]^Q(AA''efr) and log]^Q(m.5) within the displayed range. 
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Figure 6: (Top) the parameter space {ms^Us) used for comparison with Dodelson, Mel- 
chiorri & Slosar in the particular case of thermal relics. The thin bands delimited by 
blue/dot-dashed lines show regions of equal fs (assuming ujaui = 0.11 ± 0.01); the ma- 
genta/dotted lines correspond to fixed values of the velocity dispersion today; horizontal 
red/solid lines to fixed ANes- (Bottom) same with, in addition, the regions allowed at the 
68.3% (Icr), 95.4% {2a) and 99.7% (Sa) C.L. by our cosmological data set, in a Bayesian 
analysis with flat priors on log]^o("^s) ^-nd Ug within the displayed range. 
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